Leading Corrections to Finite-Size Scaling for Mixed-spin Chains 
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We identify the leading corrections to Finite-Size Scaling relations for the correlation length and 
twist order parameter of three mixed-spin quantum spin chains for the critical feature that develops 
at, 9 = TT, corresponding to a change in the topological realisation of the groundstates. 



The effective model that describes the groundstate of 
quantum spin chains is the 2d 0(3) model with a nonzero 
^?-term [ll]. Whilst it is clear that the groundstates of 
quantum spin chains can be related to the renormaliza- 
tion group fixed points of the SU(2) WZW model 0, 
it is less clear how the phase-space of just one nonper- 
turbative system can be defined from the perspective of 
renormalization. What we would like to understand, via 
the 2d 0(3) model, is how the values of the couplings in 
the action (spin stiffness, spin wave velocity, etc.) influ- 
ence the groundstate. This would tell us how the renor- 
malization group flow of the SU(2) WZW can be related 
to quantum spin chains on a physical scale. However, 
the value of these couplings is dependent on the nonlocal 
fluctuations of the system, for which we either have to 
make some form of ansatz, or determine phenomcnolog- 
ically. Conversely, a Monte Carlo calculation is a very 
good way of generating a system with nonlocal fluctua- 
tions. Locally, the interaction couplings of numerical lat- 
tice system, such as a Quantum Monte Carlo simulation, 
are well-defined. We could therefore imagine relating the 
coupling constants of the 2d 0(3) model to the lattice in- 
teraction couplings to quantify the nonlocal fluctuations. 
However, although it is meaningful to vary the value of 
the couplings arbitrarily in numerical methods, not all of 
the corresponding physics is renormalizable as a contin- 
uum theory. The difficulty is that, 6, cannot be varied 
continuously between, = 0, and, 9 = n, to describe the 
groundstate of the 0(3) model. There are nonintegrable 
singularities in the model that prevent this 

Several numerical studies have been performed previ- 
ously that identify the critical regions associated with the 
topological rearrangement of the groundstate in quantum 
spin chains. This critical behaviour occurs (in the ab- 
sence of external fields) when, 6, changes value between, 
6 = 0, and, 6 = tt, via changes in the nonperturbative 
vacuum What has not been identified in these 

previous numerical studies is whether or not the non- 
perturbative behaviour constitutes a continuous topolog- 
ical phase transition in, 9. In principle there can be 
a second order transition at zero temperature in quan- 
tum spin chains without contradiction with the impor- 
tant Mermin- Wagner theorem [§] . This theorem excludes 
the development of long-range order in finite temperature 
quantum spin chains, and therefore the development of a 



genuinely massless phase corresponding to a broken sym- 
metry, but not at zero temperature. What we are able 
to do via Quantum Monte Carlo analysis is to keep the 
volume of the lattice system fixed in physical units, but 
to vary the size of the nonlocal fluctuations by changing 
the number of lattice sites in the physical volume. The 
purpose of this new study is therefore to try and quantify 
the renormalization group flow around, 9 = ir, in terms 
of nonlocal fluctuations, to deflne the physical scale of 
the renormalization group flow. 



I. MIXED-SPIN MODEL 

In this article we investigate the nature of the ground- 
state critical feature corresponding to, 9 = n, hy study- 
ing chain systems that are amenable to numerical eval- 
uation: ones for which we can circumvent the numerical 
complex-action problem pLO| . Speciflcally, we consider 
AFM mixed chains for the Hamiltonians with spins S'* 
and S*^ arranged in cells of period-4. We choose with- 
out loss of generality Jaa = Jbh = 1, and the coupling 
anisotropy a = Jab/>/aa, L being the chain length. We 
vary the coupling interaction strength between spins to 
effectively vary, 9, through, a. 
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The periodicity in ([1]) is a necessary condition for the 
constraint equation used to define the groundstate given 
in terms of a nonlocal singlet operator [l| . A detailed 
picture of all possible realisations of these singlets for 
mixed-spin chains is given in p^ . Specifically, if three 
different couplings are used for the periodic cell then 
it is possible for these singlets to be realised between 
ncxt-to-nearest-neighbour sites, but here we restrict our- 
selves to the case of two AFM couplings and consequently 
have nearest-neighbour singlets. Our reason for looking 
at mixed-spin chains, rather than the more generalised 
bond-alternating chain p^ . is that there are differences 
in these groundstate parametrisations, as a function of 
spin magnitude, that allows us to identify important 
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nonlocal trends in the numerics. In these groundstate 
parametrisations the nonperturbative condition, 9 = n, 
arises when the singlet source term is vanishing At 
this point for a bond-alternating chain. 



(2) 



where, Sa, is the spin magnitude, and for a mixed-spin 
chain, 

— = 2SaSb{Sa + Sb)/ (s^ + + 2SaSh Jaa/^ab)- (3) 
ZTT 

In the first case, for bond-alternating chains, Jab/Jna = 
1 is always a solution for half-integer systems, and for 
integer systems the critical ratio of, Jab/^aa, tends to 
zero with increasing, Sq. In the second mixed-spin case 
the critical ratio of couplings always differs significantly 
at successive critical regions within a given model, at 
6 — TT, Stt, Stt... , etc.. The value of, 9, associated with 
critical behaviour therefore has an implicit dependence 
on the spin magnitude for mixed-spin chains, which is 
not as evident for bond-alternating chains. 

Although in these treatments the groundstate 
parametrisations can be used to identify the number of 
critical points in, 9, these parametrisations cannot be 
used to properly define local expansions in, 9, about the 
critical point at, 9 = tt, in physical units. The reason 
for this is that a simple ansatz has to be defined in order 
to be able to quantify the nonlocal fluctuations of the 
vacua. However, by definition this is a nonperturbative 
quantity, because the point, 9 = is asymptotically free 
[3|. Consequently there are major differences between 
the value of the critical point in physical units in the nu- 
merics [l|-Hl and these analytic treatments because 
the scale of the nonlocal fluctuations is essentially differ- 
ent in each case. To resolve this discrepancy we really 
need to parametrise the scale dependence of the nonlocal 
fluctuations in terms of physical units. 



II. NONPERTURBATIVE RENORMALIZATION 
GROUP FLOW 

Recently, what has been proposed is that the nonlo- 
cal fluctuations lead to an additional term in the action 
of the corresponding conformal model, in the vicinity of, 
9 = TT, [iBltlfii]- This gives rise to a double Sine-Gordon 
model description of the nonlocal fluctuations. The orig- 
inal Sine-Gordon model has been presented in a lattice 
formalism in 1^ 1^. As usual, when the topological 
Berry-phase term vanishes at, 9 = n, the action is de- 
scribed by the U(l) Gaussian theory given by the pure 
Sine-Gordon model. In terms of the lattice formalism 
the renormalization group equations of this original Sine- 
Gordon model is given by. 
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where, I = logi, and, a,, and, a,', are nonperturba- 
tive or phenomenologically determined parameter, which 
relate respectively to the couplings of the Berry-phase 
term and the kinetic terms. In the original Sine-Gordon 
model, these relations in ([4]) imply that there is a renor- 
malization of the gap in the vicinity of the critical point. 
This has be treated explicitly by via perturbative expan- 
sion in [l^. The difference for the double Sine-Gordon 
model is that the treatment has two terms correspond- 
ing the topological phase, which have different prefactors. 
Thus, in addition to the renormalization of the gap, the 
inclusion of nonlocal fluctuations will lead to the topolog- 
ical Berry-phase being effectively made anisotropic via a 
Lorentz boost. The leading correction to (|4|) should then 
be a relatively simple factor, which accounts for this ef- 
fective anisotropy. For the second order scaling in ^ 
this anisotropy will mean a logarithmic correction to the 
free energy, defined in lattice units. 

Of course, both linear and logarithmic corrections to 
FSS can arise in the pure Sine-Gordon model from con- 
formal symmetries. Quantum systems inherently have a 
conformal charge, c, that can be derived from their asso- 
ciated Virasoro algebra [2^ [2l| . For the Gaussian system 
at the renormalization group fixed point (with c = 1) this 
will lead to an additional term in the free energy with an 
explicit temperature dependence. Similar universal cor- 
rection contributions can also arise from the topology of 
the Lattice system: in it was argued that bound- 
ary conditions can lead to linear corrections to FSS from 
the conformal block. These are all potential indepen- 
dent sources for corrections to scaling, from the source of 
scaling we wish to focus on. From Haldane's conjecture 
all quantum spin chains mapping to the point, 9 = n, 
should share these universal conformal corrections and be 
of the same universality class. In this article, however, we 
arc specifically interested to identify behaviour relating 
to the extra term in [l5l|[l6l|. The contributions coming 
from this term are nonuniversal, being dependent on the 
nonlocal fluctuations of, 9. For the mixed-spin models 
we are considering, however, the critical anisotropy, a, is 
strongly dependent on spin magnitude for the different 
transition between topological sectors ([3]), even though, 
c = 1, for all of these critical points. A strong spin mag- 
nitude dependence, and logarithmic corrections, should 
therefore be indicative of the nonuniversal FSS correction 
due to the nonlocal fluctuations in our study. 



III. NUMERICAL MEASUREMENTS 

We consider the FSS behaviour of two thermodynamic 
indicators: the correlation length, evaluated via a second- 
moment estimator method [23|, and the twist order pa- 
rameter [T^]. The second-moment estimator for the cor- 
relation length, ^, is defined through the Fourier trans- 
form of a series, and potentially, this series is divergent 
away from zero temperature. Also, since we expect the 
critical behaviour to be of the SU(2) WZW universality 
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class, potential systematic errors arc indicated from ana- 
lytic results for half-integer chains at finite temperatures 
[2J]. These results identify the presence of logarithmic 
corrections to the finite-temperature scaling relations, al- 
though, this should not effect our identification of nonuni- 
versal corrections to FSS. The twist order parameter, zl, 
can be considered as a generalisation of the string order 
parameter to systems of broken translational invariance, 
and it is an exact order parameter for the groundstate 
singlet [13 • The asymptotic form of, z^, is given by, 

ZL = (vl/(n)|[/|*(")) = - 0(1/L)], (5) 

where, n, counts the number of groundstate singlets, 5*, 
of a given topology per periodic cell, and the operator, 
U, is defined by. 



U = exp 



(6) 



where, Sj, is the z-component of spin at lattice site, j, [j]. 
This order parameter changes sign as it passes through 
the critical topological region. From the asymptotic defi- 
nition in ([5]) we may expect generic linear corrections to, 
zl, that potentially affects the FSS. However, this again 
should not prevent us from identifying the presence of 
non-universal corrections to FSS. 



IV. RESULTS AND ANALYSIS 



We have investigated three mixed-spin chain systems 
with the periodic cells | — f^l^l: l^l^l^l 
i — i — 1 — 1 which we label A, B and C. Fits to the 
FSS trial functions were extracted over the numerically 
accessible temperature regimes for each spin system. Be- 
low the smallest temperatures presented in the tables we 
found that the continuous-time Quantum Monte Carlo 
method we used became inefficient to update, and above 
the largest temperatures the data became too smooth to 
fit. This is consistent with the Sine-Gordon model renor- 
malization group flow in which indicates that above 
a certain temperature scale the flow of the gap should go 
to a finite value above the gapless point. From the possi- 
ble _topological realisations of the groundstate identified 
IJI, we confirmed that model A has two critical 
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regions, whilst B and C have each one, although there 
is strong quantitative disagreement about the location of 
these transitions. We believe this is due to the renor- 
malization scale dependence of the nonlocal fluctuations, 
which we now explicitly treat in this analysis. For our 
numerical measurements we generated 100,000 indepen- 
dent conflgurations per Lattice ensemble for system sizes 
L=16, 32, 64, 128, 256, and 512. The critical value of 
the coupling anisotropy, q;c,l, was determined on finite 
chain lengths from interpolation using the Levenberg- 
Marquadt method to find the maxima of, ^, for the cor- 
relation length, and the point at which, z^ = 0, for the 
twist order parameter. 
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FIG. 1: Correlation length maxima, ^mai, versus inverse tem- 
perature, P, as a function of chain length, L. 



Our basic assumption is that the scale of the nonlo- 
cal fluctuations scales linearly with spatial length in the 
method we have used. We verify this in Fig.l, where we 
plot the length dependence of the maxima of, ^. When 
either the inverse temperature, /3, is increased at fixed 
chain length, L, or vice versa there is a departure from 
linear (second order) scaling when either scale approaches 
the magnitude of the correlation length. However, with 
fixed ratio between the two scales the correlation length 
can always be made to scale linearly with, L. Therefore 
there is a direct correlation between the chain length and 
the size of the system in physical units of, /3. 

We extracted the shift exponent. A, via fits to. 
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a/L + b/L^ 



(7) 



where, ac, a, b, and. A, are fitted constants. The results 
are presented in Tables 1 and 2. To summarise: no stable 
fit was obtained for the correlation length measurement 
of the second transition point in A which corresponds 
to, 9 = Stt, and the exponents all appear highly nonuni- 
versal. We anticipate both universal linear and nonuni- 
versal logarithmic corrections to scaling due to, respec- 
tively, the conformal symmetries and the nonlocal fluctu- 
ations in our analysis. Although we can generate lattice 
ensembles to correspond to a range of different nonlo- 
cal fluctuation scales, it is difficult for us then to make 
a direct comparison between nonuniversal and universal 
corrections. The problem is that in the continuous-time 
method varying the scale of the nonlocal fluctuations also 
changes the renormalization scale of the gap, from 
because the size of the lattice is defined in physical units, 
rather than mmrber of lattice sites. Although the size 
of the conformal correction should be fixed, the lattice 
volume is fluctuating. The way we treated this was fol- 
lowing [2^ , to look at the relative strengths of universal 
and nonuniversal corrections. Since the lattice units are 
nonuniversal, we instead fix the nonuniversal exponent. 
The nonuniversal exponents, from Tables 1 and 2, range 
between 1.4-3.1, and so we fix the value of the exponent 
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0.01 




0.48224(15) 


1.23(24) 




zl 


0.48360(48) 


1.32(5) 




ZL 


1.3126(14) 


1.41(9) 


0.025 




0.48279(21) 


1.35(32) 




ZL 
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1.3146(7) 


1.47(6) 


0.04 




0.48314(18) 


1.51(17) 




Zl 


0.48189(47) 


1.25(4) 
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1.3146(5) 


1.38(6) 



TABLE I: Model A, 
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0.01 


Zl 


0.62132(5) 
0.62107(49) 


2.32(47) 
1.59(5) 


0.01 


0.76244(5) 
0.76246(17) 


2.72(14) 
1.7(1) 


0.025 


ZL 


0.62125(17) 
0.62141(17) 


1.62(47) 
1.57(4) 


0.025 


0.76289(8) 
0.76389(17) 


3.1(3) 
1.55(6) 


0.05 


ZL 


0.62184(33) 
0.62256(52) 


1.32(29) 
1.44(4) 


0.04 


0.76440(13) 
0.76508(29) 


2.78(38) 
1.58(7) 


0.075 


ZL 


0.62211(65) 


1.48(5) 


0.05 


0.76526(21) 


2.41(33) 



TABLE IL Models B, |-|-i-i,andC, |-|-1-1 



to be, A = 2, to make the comparison, via the constrained 
fit, 
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FIG. 2: Linear correction term, a, and logarithmic correction 
(8) term, 6, for the correlation length, ^, versus temperature, T. 



Tlie quahtative comparison of the hnear, a, and iog- 
arithmic corrections, b, is presented for the correlation 
length of models B and C in Fig. 2. From our fits, a, is 
significantly different from zero for B and C, and in addi- 
tion there is an unexpected temperature dependence for 
B. Similar behaviour is seen for, 6, although the size of 
correction is considerably larger, supporting our premise 
of a leading logarithmic correction. Since the tempera- 
ture is fixed in physical units for each measurement, this 
unforeseen temperature dependence is most likely due to 
the systematic errors in the improved estimator method 
for the correlation length, since it is largely absent for the 
corresponding Twist order parameter plots in Fig. 3. The 
FSS corrections are strongly dependent on the spin mag- 
nitudes for the specific models. A, B and C, which also 
suggests that the leading corrections to scaling at the, 
c = 1, point are nonuniversal. The Twist order parame- 
ter measurements of Fig. 3 are generally more stable than 
the corresponding correlation length measurements, and 
we are able to more distinctly see the large separation in 
values between the size of the linear and logarithmic cor- 
rections. The twist order parameter is perhaps a more 
consistent with our expectations, because it is a more 
direct measure of the singlet state than the correlation 
length second-moment improved estimator. 



V. SUMMARY 

In this article we have identified the leading correc- 
tions to FSS of three mixed-spin chains, via Quantum 
Monte Carlo analysis, in the vicinity of the critical point 
at, 9 = TT. The scaling picture presented is one defined 
in terms of, /?, and the effective value of, 6, determined 
from the fluctuations of the lattice units, and the transla- 
tionally varying interaction coupling, Jab- We have iden- 
tified that the leading corrections are logarithmic and 
nonuniversal, dependent on both the spin magnitude of 
the system and the nonlocal fluctuations of the lattice 
ensemble. Measurements have been compared between 
a second-moment estimator measurement for the corre- 
lation length, and the twist order parameter, and it was 
found that the latter were the more numerically stable. 
Our numerics conflrm the effective renormalization re- 
gion picture developed in [l^ [l^ of the addition of non- 
local fluctuations to the scaling picture of the c = 1 con- 
formal treatment. However, it is not clear that we can 
determine the dependencies on, 0, and, /3, from our mea- 
surements sufficiently to be able to define the correspond- 
ing lattice /3-function. The reason for this is that the 
continuous-time method we have used is defined in fixed 



5 



1 

0.5 


-0.5 
-1 

25 
20 
15 
10 

5 



-5 - 
-10 



"T 1 r 



"T 1 — 

A — X- 
B B- 

C - -e> 



_l I L_ 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 
T 



"T 1 1 1 1 r 



A -x- 
B B- 

C - -e> 



_l I I I I L_ 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 
T 



physical units, but with a fluctuating space-Euchdcan 
time volume. Consequently, the lattice units of our FSS 
are phenomenological, and it is difficult to obtain an un- 
renormalized value for the central charge. It is therefore 
difficult for us to make a comparison of the unrenormal- 
ized nonuniversal and universal corrections of the system, 
since the lattice units are themselves nonuniversal. How- 
ever, we have been able to make a qualitative comparison 
of the size of the universal and nonuniversal correction 
terms by picking an arbitrary common renormalization 
scale for the scaling functions. Our numerical results for 
the nonuniversality indicate that the contribution from 
different topological sectors, defined in terms of the topo- 
logically distinct groundstate realisations, is crucial to 
understanding the renormalization group flow of numer- 
ical Quantum Monte Carlo data. Although our specific 
analysis has been performed for quantum spin chains in 
the vicinity of the zero temperature fixed point of the 
double Sine-Gordon model, we anticipate our findings 
may be of general applicability to low-dimensional quan- 
tum spin systems, given the ubiquitous character of the 
0(3) model in describing topological quantum effects. 



FIG. 3: Linear correction term, a, and logarithmic correc- 
tion term for, b, for the Twist order parameter, zl, versus 
temperature, T. 
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